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Abstract. - We propose two lattice models in one dimension, with stochastically hopping particles 
which aggregate on contact. The hops are guided by "velocity rates" which themselves evolve 
according to the rules of ballistic aggregation as in a sticky gas in continuum. Our lattice models 
have both velocity and density fields and an appropriate real time evolution, such that they can 
be compared directly with event driven molecular dynamics (MD) results for the sticky gas. We 
demonstrate numerically that the long time and large distance behavior of the lattice models 
are identical to that of the MD, and some exact results known for the sticky gas. In particular, 
the exactly predicted form of the non-Gaussian tail of the velocity distribution function is clearly 
exhibited. This correspondence of the lattice models and the sticky gas in continuum is nontrivial, 
as the latter has a deterministic dynamics with local kinematic constraint, in contrast with the 
former; yet the spatial velocity profiles (with shocks) of the lattice models and the MD have 
striking match. 
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, Introduction. — A gas of ballistically moving parti- 
cles which suffer completely inelastic collisions amongst 
themselves and thereby aggregate, is weW known as the 
^ Sticky gas. The motion of the particles are deterministic 
g,nd the collisions follow the momentum and mass con- 
servation laws. The problem becomes statistical as the 
initial velocities of the particles are assumed to be ran- 
dom. The system exhibits large scale density clustering 
^^"with a growing coarsening length scale, and appearance 
T— (of shocks in spatial velocity profile. The connection be- 
tween this problem and shock dynamics of a Burgers fluid 

* * 

^ in the high Reynolds number limit was established in [1] . 
? The sticky gas problem has some relevance for interstellar 
rS matter, as similarities between self gravitating matter and 
j^ijaertial Burgers fluid was shown in [2]. In the context of 
terrestrial dissipativc granular matter, a mean field anal- 
ysis of the model with several scaling predictions was in- 
troduced in 13]. The model was exactly solved for several 
statistical functions in [3J[S]. 

Apriori the sticky gas, which is a gas of particles with 
zero restitution coefficient (r) may appear to be rather ar- 
tificial, since real granular gases all have < r < 1. In 
fact in reality, r is also dependent on the relative velocity 
Wroi of coUisional impact — for the Wrci — > 0, r H> 1 [SJ[7]. 
Very interestingly it has been shown that for intermediate 



times (but not very large [5]) such a realistic granular gas 
with any r (0 < r < 1) behaves statistically similarly as 
a sticky gas (r = 0) [Hllin]- This remarkable universality, 
hitherto is only supported by molecular dynamics (MD) 
and not proved analytically. In particular, since the map- 
ping between the particle-dynamics of the sticky gas and 
shock-dynamics of the Burgers fluid is known [T] and sev- 
eral properties of the sticky gas and the Burgers equation 
has been derived analytically [4l[5j[Tl], one hopes that the 
speculation of the inviscid Burgers equation being the cor- 
rect continuum limit of the finite r granular gas should be 
analytically provable. 

In certain non-equilibrium problems like diffusion lim- 
ited reactions [T^HT^ , a quantum field theoretic formalism 
has been develop to systematically connect the stochastic 
dynamics of appropriate lattice models to their respective 
continuum Langevin dynamical limits. Motivated by the 
latter works we wonder if a similar method can be fol- 
lowed for the granular gas problem discussed above. Can 
one invent suitable stochastic lattice models which may 
represent the microscopic ballistic granular gas problem? 
The proposal is rather ambitious for r 7^ 0, but one may 
start with the simpler sticky gas (i.e. the r = case). In 
this paper, we propose two stochastic lattice models (with 
single and multi-particle occupancies per site respectively) 



p-1 



Supravat Dey et al. 



which are shown to have the same statistical properties, at 
long times and large distances, as the sticky gas. Statisti- 
cal distribution functions of the mass and velocity, spatial 
density-density correlation functions, and velocity shock 
profiles, for our lattice models are shown to have strik- 
ing match with exact results and MD simulation results 
of the sticky gas. Of the two models we discuss below, 
the single particle occupancy model would seem a natu- 
ral counterpart of the sticky gas. On the other hand, the 
second model with multiple site occupancy and a delayed 
aggregation would apriori seem a bit different — although 
in the limit of a tuning parameter we show that it reduces 
to the single occupancy model. Explicit study of a multi- 
ple particle occupancy model is motivated by the fact that 
the field theoretic literature [T2j[l3] have generally dealt 
with multiple occupancy models. Furthermore they are 
also natural for future extension of our study to the r ^ 
case, where particles do not stick after collisions. Thus 
our choice of the models leave open the scope of coarse 
graining in future, and obtaining the desired continuum 
limit. 

Although lattice models have been studied earlier [T51 - 
119] for the freely cooling granular gas, to understand the 
velocity ordering and formation of spatial shocks they had 
many limitations. Firstly, these earlier models only had 
lattice velocities which evolved stochastically and they in- 
volved no actual particles hopping. As a result although 
key insights on the behavior of the velocity field could be 
obtained, the "density field" could not be tracked. More- 
over the time used for characterizing evolution of the sys- 
tem was "collision number" , and matching to real time 
evolving MD involved arbitrariness. Finally kinematic 
rules had to be specified separately on the lattice mod- 
els to avoid collisions disallowed by the ballistic dynamics. 
In the models that we present in this paper, these short- 
comings have been overcome. We have actual particles 
hopping and so both the velocity and mass density fields 
can be tracked. The time evolution is done using prescrip- 
tions of exact stochastic simulation of Master equation, 
such that it compares directly with time in MD simula- 
tions. Interestingly the stochastic dynamics in our models 
naturally take care of one type of kinematically forbidden 
collisions, but violates another type (a slower particle can 
catch up with a faster one from behind). In spite of this 
lack of strict adherence to kinematic constraints we show 
that the results of velocity shock profiles are neverthe- 
less almost identical to MD for identical initial conditions. 
Thus we claim our lattice models to be minimally complex 
models necessary to represent the sticky gas system, and 
may hope (as discussed above) that in future they may be 
appropriately coarse-grained to derive the inviscid Burgers 
equation. 

Below we begin by first defining the single and multiple 
occupancy lattice models SOSG and MOSG, respectively. 

Model I: Single occupancy sticky gas (SOSG). 

The system has massive particles on a one-dimensional 



periodic lattice of size L. An occupied site i at time t 
(where i can take values 1,2, L) has two non-zero vari- 
ables: mass m.i and "velocity" Vi. The magnitude of the 
velocity Vi specifies the rate of hopping of the particle to 
its nearest neighbors, and its sign decides the direction of 
hopping. For an unoccupied site we designate = and 
Vi = 0. Wc initialize the system by having nii = 1 Vi and 
assigning random velocities Vi drawn from an uniform box 
distribution over the range ( — 1, +1). The system evolves 
in time via hopping and instantaneous aggregation events. 
Let the sites i and its neighbor j have variables {mi,Vi) 
and {mj,Vj) respectively at time t. The particle of the 
site i may hop to either of its neighbors j(= i + l,i — 1) 
depending on sgn(vi){= -1-1,— 1), and with the rate \vi\. 
If at time t + 6t the particle at the z*'' actually hops to 
the j then the masses and the velocities change in the 
following way: — >■ 0, — >■ 0, rrij rrii + rrij, and 
Vj — ;> {rriiVi + rrijVj) / {nii 4- nij). The latter update of Vj 
follows the momentum conservation rule of particle colli- 
sion with zero restitution coefficient. 

To implement the above dynamics, in our stochastic 
simulation, we choose an event of hopping of a particle 
at site i with probability |wi|/r [5DH53], where F = \ vi\ 
is the sum of all hopping rates possible in the system at 
the time t. The enactment of the event is associated with 
a time increment 5t, which we choose to be the reciprocal 
of the total rate F, i.e 5t = 1/F. As F is time dependent 
in this problem, 5t is also time dependent. We note that 
this choice of 5t is the average value of the truly random 6t 
drawn from its distribution function Fexp(— F(5t) |20H23) . 
Although using a random 5t would have been more appro- 
priate in principle, for very small rate F as we have at large 
times, sampling its full distribution becomes very tedious 
numerically. Using instead the average value 5t = 1/F 
saves numerical effort considerably, without compromis- 
ing accuracy (as we have checked). 

Model II: Multiple occupancy sticky gas 
(MOSG) . — In this model we also have a system of mas- 
sive particles on an one-dimensional periodic lattice of L 
sites, except that now multiple particles can occupy any 
site. For a site i occupied by a*'' particle, we associate two 
variables mi^a and Vi^a- If the number of particles at site i 
is TT-i, then a can take values 1, 2, rii. As in SOSG 
denotes the rate of hopping while sgn{vi^a){^ +1, —1) de- 
cides to which neighboring site j = {i + l,i — 1) the par- 
ticle may hop. The system evolves via hopping followed 
by aggregation, except that in MOSG unlike SOSG, ag- 
gregation events are not instantaneous. At any site i, a 
pair of particles a and /3 may aggregate with a rate A. 
In the limit of A — !> oo, aggregation events becomes in- 
stantaneous leading to maximum occupancy of the sites 
being 1 — thus MOSG maps back to SOSG. We initial- 
ize the system by having n,; = 1 and m^.Q ~ 1 Vi, and 
assigning random velocities Vi^a drawn from an uniform 
box distribution over the range (— 1,+1). After aggrega- 
tion of two particles a and /3 at site i, nii^a nii^a + Tni^p, 
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Vi,a {mi^aVi^a+rrii^pVi^p) / [nii^a+mi^jj), and n,; n^-l. 

As in SOSG the stochastic simulation is implemented 
by choosing events of hopping of particles and their ag- 
gregation randomly with probabilities \vi^a\/^ and A/F 
respectively. Here T = J2i,a ki.^l + Si " 1)^/2 is 
the sum of all hopping and aggregation rates possible in 
the system at the time t. Note that the factor ni{ni — 1)/2 
comes from the number of possible collisions among any 
pair of particles at site i. The increment of time between 
two events is chosen to be St = l/F. 

As noted above, both the models SOSG and MOSG are 
more realistic in having on one hand both local mass den- 
sity and velocity fields, and on the other hand realistic 
time increments St. In both these respects we overcome 
serious drawbacks of the earlier lattice models [TSUTTlfTQ] 
of dissipativc gases in making contact with molecular dy- 
namics and continuum theories. At the same time the 
models being stochastic and not having strict kinematic 
constraint (as will be discussed below), leave the curi- 
ous question open as to whether their long time and large 
wavelength properties would match with that of the de- 
terministically evolving ballistic aggregation model. We 
now proceed to show numerically that such a non-trivial 
correspondence is indeed exhibited by both the models. 

Results. - For both SOSG and MOSG lattice mod- 
els we use system sizes L = 50000. As mentioned before 
the hopping rates Vi (in SOSG) and Vi^a (in MOSG) are 
initially chosen to be random and drawn from the uniform 
box distribution over range (-1,1). In MOSG we have used 
the aggregation rate A = 5.0. The unit of time t in our 
lattice simulations is inverse of the unit of rates — as rates 
are chosen to be dimensionless, the time is also dimension- 
less. Similarly the space x in units of the lattice spacing 
is taken as dimensionless. 

The choice of the value of A = 5.0 is based on the ob- 
servation of energy decay in MOSG for various values of 
A shown in Fig. [TJ The scaled energy e(t) for the MOSG 
model with different A, tend to that of the SOSG model 
for increasing A. For large t both the models follow the 
classic energy decay e{t) - t^^/a ^ell known for 

the ballistic sticky gas. In particular for A = 5.0 we can 
see from Fig. [T] that for times t > 10^ the MOSG and 
SOSG curves are completely merge, implying that we are 
safely in the universal sticky gas regime. Interestingly the 
irrelevance of A and the universality of e(t) at large t, is 
reminiscent of a similar irrelevance of the values of restitu- 
tion coefficient in connection to energy decay of granular 
gases in one-dimension j9j . 

Exact analytical results in terms of explicit functions 
and integrals for the mass and velocity distribution func- 
tions of the sticky gas are available [4l[5]. On the other 
hand the mass density-density spatial correlation function, 
although not available as an explicit analytic form, can 
easily be obtained from event driven molecular dynamics 
(MD) simulation of the sticky gas [T^. Below we first pro- 
ceed to match the numerically obtained mass and velocity 



distributions from the SOSG and MOSG models with the 
exact analytical results. Then we compare the mass den- 
sity correlator of SOSG and MOSG models with that of 
MD. Finally, even the microscopic spatial velocity shock 
profile of the two lattice models are shown to be almost 
identical to that of the MD. 




Fig. 1: Log-log plot of Energy e{t) = E{t)/Eo with time t, 
where E{t) is the total energy at time t and Eq = E{t = 0). 
Data points with lines joining them are for the MOSG model 
with different values of A indicated as labels. The solid line 
is for the SOSG model. The dashed line has slope -2/3 and 
serves as a guide to eye. 

Comparison of the lattice model results with the analyt- 
ical results. The exact analytical solution of the sticky 
gas problem [?1[S] relies on reducing the calculation of the 
distribution functions to a summing of Brownian paths 
with parabolic constrains in momentum space. The dif- 
fusion constant D associated with the effective Brownian 
motion is exactly equal to the inverse of the variance f3 of 
the initial momentum distribution of the particles in the 
sticky gas [5] . The following distribution functions for the 
scaled mass M = m/t^/'^ and scaled velocity V = vt^/'^ 
were derived in |5] as follows: 



Mi(M) = pMI{M)H{M) 
JiiiV) = 2^0 ^ dMMI{M)x 



(1) 



J{V - M/2)J{-V - M/2) (2) 



where, 
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(3) 



In Eq. [2]Ai is the Airy function [H] and -uJk (k=l,2,3,...) 
are its zeros. 
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The only parameter appearing Eqs. [TJ [51 and [3] is /? 
which has to be fixed to match these formulas with our 
simulation results. In both the lattice models the initial 
mass of every particle is unity, and hence the variance of 
the momentum distribution is same as that of the velocity 
distribution. The variance of the uniform box distribution 
of velocities over the range (-1,1) is /3 = 1/3. Using the 
latter value of f3, we evaluate the integrals H and using 
the "quadgk" function with absolute tolerance — 10~^^ 
in MATLAB. The evaluation of the sum T{M) was done 
using our own code — for large M (M > 0.1) we kept up 
to 6000 terms in the sum, while for small M [M < 0.1) 
we used the asymptotic formula I{M) ~ 1/ y^2T:/3M^ [S]. 




Fig. 2: Log-log plot of P{m,t)t'^^^ versus M = m/t'^''^ for two 
different models at two different times. Circular symbols (O, 
• ) represent the SOSG data, while triangular symbols (A, 
a) represent the MOSG data. Empty symbols are data for 
t = 8 X 10^ and solid symbols are data for f = 16 x 10^. The 
solid line is the exact analytical formula of Eq. [T] 

The mass distribution function P{m, t) is related to the 
scaling function (M) (Eq. [T]) as follows: 

P{m,t) = ^fi,{M). (4) 

In the simulation of the SOSG model at any given time t, 
we drew the frequencies of the non-zero integer masses 
sweeping over all the sites of the lattice. In contrast, 
in the MOSG model the frequencies of the total mass 
(SaLi iTT-i.a) of every occupied site were drawn, again 
sweeping over the lattice sites. Normalizing the frequen- 
cies by L and averaging these further over several random 
initial conditions (^ 10*) we finally obtained the mass dis- 
tribution function P{m,t). The latter distribution func- 
tions obtained for the stochastic SOSG and MOSG models 
at two different times namely t = 8x 10^ and t = 16x 10^ 
are scaled (Eq. |4]) and plotted in Fig. [2l The exact scal- 
ing function /ii(M) enumerated by the method described 
above is plotted against the data of the lattice models. 
The spectacular match of the Monte Carlo data and the 
exact formula, without any parameter fitting, gives the 



first evidence to our claim that the lattice models SOSG 
and MOSG are correct representatives of the continuum 
behavior of the sticky gas. 

We now proceed to look at the velocity distributions 
for the lattice models. The normalized distribution of our 
interest scales as follows: 

Q{v,t)^t'^^h{vt'^^), (5) 

where fi{V) = fli{V)/ p-i{V), and iIi{V) is given by 
Eq. [2] Unlike the integer masses m, the velocities are real 
numbers, so a numerical binning is required to obtained 
a distribution function. For SOSG, at any time t wc find 
the frequencies of velocities of occupied sites falling within 
coarse-grained bin widths of magnitude 0.001. For MOSG, 
the frequencies of average velocities (X^qLi Wi,a/'^i) of oc- 
cupied sites falling within different bins (each of width 
= 0.001) were obtained. We normalize these frequencies 
by the total number of occupied sites and the bin width. 
Since the tail of the velocity distribution is of great in- 
terest, to obtain high precision, we had to average over 
relatively much larger set of initial conditions in this case 
(namely ~ 10^). 















10-^ 










Q(v, t) 
















1.9 


2.3 y 


2.7 ^\ - 



-2 -1 y 1 2 



Fig. 3: Log-linear plot of Q{v,t)/t^^^ versus V = vt^^^ for 
two different models at two different times. Circular symbols 
(O, • ) represent the SOSG data, while triangular symbols 
(A, a) represent the MOSG data. Empty symbols are data 
for t = 8 X 10^ and solid symbols are data for t = 16 x 10^. 
The solid line is the fi{V) using exact analytical formula for 
/ii(F)(Eq. The dotted line is a Gaussian curve with zero 
mean and variance = 0.26 which fits only the small V data. 
Inset: A zoom into the tail region of Q{v,t) /t^^^ versus V (for 
positive V) to show clearly the deviation of the data and the 
exact formula from the Gaussian form. 

The velocity distribution data obtained from the lattice 
models SOSG and MOSG at two different times namely 
t = 8 x 10^ and t = 16 x 10^ are scaled (Eq. [5]) and 
plotted in Fig. [3l The exact curve of fi{V) (obtained 
using p,i{V) of Eq. ^ shown with a solid line in Fig. [3] 
passes all the way through the Monte Carlo data. Again 
the spectacular match of the data and the exact result 
without any parameter fiddling is to be noted. A Gaussian 
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curve shown in dotted line fits the data for small V, but 
the tails of the data depart from it. 

The asymptotic form of the Eq. [5] for large V, i.e. the 
tail, is not known exactly. In [S] the tail was shown to be 
bounded as /2i (F -> oo) < constant exp(-/3|y|3/6 - 
{l3/2)'^^^\V\uji). The inset of Fig. [3] shows the Monte 
Carlo data near the tail; the deviation from Gaussianity 
and match with the exact Eq. [5] are conclusive. How- 
ever, to determine the precise asymptotic decay, we will 
require data for much larger range of V. At this point we 
would like to recall that finding no n- Gaussianity in nu- 
merical study of freely cooling granular gases has been an 
extremely challenging task in earlier works. For example 
in one-dimensional granular gas with finite restitution co- 
efficient [5], although the system was argued to approach 
the sticky gas limit asymptotically, the numerical data for 
the velocity distribution did not clearly show deviation 
from Gaussianity. We have checked that the computation 
time involved in MD is very huge to attain the neces- 
sary accuracy to demonstrate the deviation of the tail of 
the distribution from the Gaussian form. Interestingly in 
a very different context of quasi-elastic limit of granular 
gas, 'stationary' velocity distribution with tail of the form 
exp(— c|l/|^) have been numerically demonstrated |25| but 
that does not directly address the concerned limit (r = 0) 
of this paper. Given this background it is significant that 
wc have succeeded in demonstrating the deviation from 
Gaussianity of the tail of the velocity distribution with 
comparative ease in the Monte Carlo simulations of the 
lattice models SOSG and MOSG. 

Comparison of the lattice model results with MD results. 
Event driven MD of ballistically moving sticky particles in 
continuum is expected to yield numerically exact results, if 
done with sufficient accuracy. The details of such simula- 
tion can be found in a recent publication |10j . Here we use 
event driven MD to study the spatial mass density-density 
correlation function and spatial velocity profile, compare 
these with that obtained from Monte-Carlo simulation of 
the SOSG and MOSG models. 

The local mass density function is defined in the fol- 
lowing fashion. For the SOSG and MOSG models a lo- 
cal site mass density rrii is defined as the total mass at 
site i. For the MD simulation the whole continuum space 
of length L is divided into L boxes and the sum total 
of masses in the z*'' box denotes the mass density rrii 
of the box [TU]. The density-density correlation function 
Cmm{x,t) = {mi{t)mi+x(t)). Here, the angular bracket 
(.) represents average over both space (i) and initial ve- 
locities, in all the three systems. In MD simulations with 
L particles, the system is initialized to have random parti- 
cle velocities drawn from an uniform box distribution over 
( — 1, 1), and unit particle masses kept at unit separations. 
In MD too we use L — 50000. Thus the initial conditions 
are same for MD and the two lattice models, although 
their subsequent dynamics are very different. 

The mass density-density correlation functions are 
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Fig. 4: Plot of Cmm{x,t) versus x for the MD simulations and 
the two different lattice models, at two different times. Empty 
symbols are data for t — %x ICp and solid symbols are data for 
t = 16 X 10^. Circular symbols (O, • ) represent the SOSG 
data, triangular symbols (A, A) represent the MOSG data, 
and the square symbols (□, ■) represent the MD data. Inset 
shows plot of Cmm.{x.t) versus x/ C{t) for the six curves in the 
main figure — the data collapse according to Eq. (6] 

shown in Fig. |4]at two different times. At both the time 
instants the data of the lattice models and the MD sim- 
ulations fall perfectly on top of each other, without any 
parameter fixing. The profile of the correlation function 
is characteristic of the sticky gas and has been discussed 
in |10j . As the system coarsens in time, Cmm{x^t) is ex- 
pected to obey the scaling hypothesis [25] when space is 
scaled with the coarsening length C{t) ^ t^^^: 



{cit) 



(6) 



As shown in the inset of the Fig. |4l there is perfect data 
collapse following the above equation. 

Unlike spatial correlation functions which capture or- 
dering and structure formation at the macroscopic level, 
microscopic information is captured by the spatial veloc- 
ity profile of a sticky gas. The velocity profile v{x) exhibit 
shocks joined by approximate linear curves [l]|9l[T9]. In 
Fig. [5] we show particle-velocity v{x) versus x for the 
MD simulation with visible shocks (where velocity jumps 
suddenly from positive to negative) as expected. It is re- 
markable that for the SOSG and the MOSG lattice models 
the velocity profile of occupied site at the same time t, are 
almost identical as the MD (Fig. provided their initial 
conditions are identical. Although in an earlier work such 
a correspondence was demonstrated |19| , there the shocks 
in the velocity profile had an inverse correspondence to the 
shocks in MD — lattice shocks were located at positions 
where MD profile was gradual and vice versa; this was a 
consequence of the fact that there were no actual motion 
of particles along the lattice. Contrary to that since we 
have actual particles hopping on lattice, we have natural 
and direct correspondence of shock profiles with the MD 
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exact analytical and MD simulation results. The lattice 
dynamics violates the strict conditions of kinematic con- 
straint, but only mildly, so that there is no observable 
effect over slightly coarse-grained spatial scale (see Fig. 

ED. 

Based on the results that we have found, a curious possi- 
bility has been opened up to derive the continuum behav- 
ior of the granular gases in future. We hope that the two 
lattice models studied here have all the necessary ingredi- 
ents sufficient to capture the correct continuum behavior 
of the sticky gas; the models should lead to the inviscid 
Burgers equation if they can be successfully coarse grained 
by a suitable technique. 



Fig. 5: Velocity profiles v{x) plotted against space x (which is 
continuous for MD and assumes integer values for the SOSG 
and MOSG models). The profiles of v{x) at time t — 10'* 
look almost identical in the three cases, for the same initial 
realization of random velocities. 

(Fig. [5]). Moreover it was emphasized in [TH] that local 
kinematic constraints had to be separately imposed on the 
lattice dynamics to ensure that collisions disallowed by the 
ballistic dynamics do not occur, and this was argued to be 
a key factor for shock formation. The kinematic constraint 
states that ii — vl > (where vl and vr denote ve- 
locities of a pair of left and right particles, respectively) 
collisions won't occur. The collisions disallowed by the 
latter rule can be classified into two types: (a) vu > 
and tjl < 0, (b) Vfi > with same sign for both. In our 
simulation since hopping in the opposite direction of rate 
V is disallowed, case (a) is automatically taken care of — 
this we would claim is indeed the key condition for shock 
formation. On the other hand, due to the stochastic na- 
ture of hopping in our models, constraint (b) is sometimes 
violated (a slower particle may catch up with a faster one) 
— yet the latter violation has no discernible effect (Fig. 
[5]) on the microscopic shock profile. Hence we conclude 
that constraint (b) is not an essential constraint. 

Discussion. — Wc have shown that two lattice mod- 
els with single and multiple particles site occupancy and 
stochastic dynamics, can exactly reproduce the long time 
and large space behavior of the ballistic sticky gas system. 
In particular, the mass and velocity distribution functions 
are identical to the exact analytically known distributions 
of the sticky gas. The density-density spatial correlation 
functions and local microscopic velocity shock profiles of 
the lattice models and MD simulations also have a perfect 
match. 

As we have noted that the current lattice models dif- 
fer from earlier lattice models in having a genuine density 
field besides the velocity field. The time update is done 
following a prescription suited for real continuous time 
evolution, allowing us to directly compare our results with 
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